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In order for sources of coherent high brightness and intensity THz and X-Ray radiation to be 
accepted by university or industrial R&D laboratories, truly compact, high current and efficient 
particle accelerators are required. The demand for compactness and efficiency can be satisfied by 
superconducting RF energy recovery linear accelerators (SRF ERL) allowing effectively minimising 
the footprint and maximising the efficiency of the system. However such set-ups are affected by 
regenerative beam-break up (BBU) instabilities which limit the beam current and may terminate 
the beam transport as well as energy recuperation. In this paper we suggest and discuss a SRF 
ERL with asymmetric configuration of resonantly coupled accelerating and decelerating cavities. In 
this type of SRF ERL an electron bunch is passing through accelerating and decelerating cavities 
once and, as we show in this case, the regenerative BBU instability can be minimised allowing high 
currents to be achieved. We study the BBU start current in such an asymmetric ERL via analytical 
and numerical models and discuss the properties of such a system. 


I. INTRODUCTION 

The next generation light sources are to be compact, 
highly efficient, have high repetition rates and high- 
brilliance radiation pulses. One of the candidates to sat¬ 
isfy all these requirements are light sources based on en¬ 
ergy recovery linac (ERLs) [ij driven by photo-injection 
sources. Linacs parameters such as: emmitance, repeti¬ 
tion rate and bunch charge, in this case are driven by 
an electron photo-injector and laser technologies. Both 
these technologies have improved dramatically in the last 
twenty years and such linac drivers capable of generating 
femtosecond pulses can be routinely bought from spe¬ 
cialised companies. To generate a high-power, high bril¬ 
liance beam either in THz or X-ray ranges, a high charge 
electron beam is required and new developments are now 
bringing Ampere class injectors to reality 0. The in¬ 
crease of the bunch charge will lead to an increase of 
photon yield and brilliance during either x-ray Compton 
scattering Q or generation of THz radiation [^. The 
power ranges of a RF power supply required to drive 
100mA (a typical current in such accelerators) are such 
that energy recovery is required to meet the demand for 
energy efficient systems. However, adding an energy re¬ 
covery stage, while increasing the beam charge and rep¬ 
etition rate, leads to the appearance of so called beam 
break-up (BBU) instabilities These instabilities re¬ 
sult in beam trajectory shifts, energy recovery degra¬ 
dation, and ultimately, termination of the beam trans¬ 
portation. The regenerative BBU instability is especially 
damaging to ERL systems and originates from parasitic 
excitation of transverse higher order modes (HOM^ in¬ 
side the cavities. The use of the same cavity @ or 
strongly coupled cavities 0,i means that the positive 
feedback between the transverse momentum imparted by 
the HOM, and hence beam displacement and the HOMs 
amplitude is readily established. A circulating beam 


through such a system results in a growth in the beam 
displacement and dephasing for each bunch. 

In this paper we discuss a single turn SRF ERL sys¬ 
tem [9I-HIII where the beam is transported through the 
accelerating section, interaction point (IP) and decelera¬ 
tion section only once. This resembles a similar recircu¬ 
lation layout as used in nuclear medicine in the form of 
the reflexotron developed at Chalk River Nuclear Lab¬ 
oratories [T^ . but has important differences as we will 
describe below. In this model, the beam is accelerated 
inside the acceleration section while in the deceleration 
section most of the beam energy is extracted and guided 
through a resonantly coupled section back into the ac¬ 
celeration section. In FIG. [H we show a schematic illus¬ 
tration of a compact source of coherent radiation driven 
by such a single turn SRF ERL. Both sections consist of 
the same number of cells but adjusted in such a way that 
insures that only the operating mode of both sections are 
fully overlapping each other creating a single operating 
mode of the cavity, while the HOMs are separated in the 
frequency domain. 

The two sections are linked by a resonant coupler, and 
here the resonant coupling means that the two sections 
are only strongly coupled at the set of frequencies which 
are overlapping eigen-frequencies of all the three com¬ 
ponents of the system, i.e. the coupling cell and both 
sections. There is still some field leakage from one sec¬ 
tion to another but as it will be shown, the effect is rel¬ 
atively small. Due to the approach described, we re¬ 
duce the possibility for the multi-pass-regenerative BBU 
feedback mechanism (where the HOM in the accelerating 
pass is driven by the field generated by the beam in the 
decelerating pass) to be established. A potential issue 
which may be caused by the asymmetry of the structure 
is that HOMs may have different electrical centres and 
the beam cannot be at the field nulls for all the HOMs. 
This theoretically may allow an on-axis bunch to excite a 
multi-pole HOM which can deflect the following bunch. 
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FIG. 1. A schematic of possible single turn ERL system. 


The paper is organised in the following manner: in sec¬ 
tion |TT] we present the model and basic description of the 
asymmetric cavity ERL, then in section Hi Al we discuss 
the model of bunch trajectories in such ERL and discuss 
the BBU instabilities. Section imi deals with analysis of 
HOMs and their start current using the RLC circuit ap¬ 
proach. Section m is dedicated to numerical modelling 
and to an estimation of BBU start currents for a pro¬ 
posed source. In conclusion (section we discuss the 
advantages of the system presented. 


II. MODEL AND BASIC DESCRIPTION OF 
ASYMMETRIC CAVITY FOR ERL 

To increase the efficiency of linac based sources of co¬ 
herent radiation, electron beam energy recovery is re¬ 
quired. Such devices are known as energy recovery linacs 
(ERLs) and they are assumed to be attractive drivers for 
compact, energy efficient sources of THz and X-ray ra¬ 
diation. One of the major issues in conventional ERLs 


is the so called ‘beam break-up instability’ (BBU) which 
significantly limits the current transported through the 
system. In ERLs, the dominant instability is the multi¬ 
pass regenerative BBU instability where a transverse kick 
to the beam is given by a higher order mode (HOM). In 
this paper we suggest a design of a single turn ERL con¬ 
sisting of an accelerating and decelerating section, which 
are resonantly coupled. 

The ERL under consideration has two axes (FIG [U . 
An electron bunch will propagate along the first axis and 
be accelerated. An electron bunch propagating along the 
second axis will be decelerated and feed energy back into 
the ERL. The cells’ shapes on each axis are tuned to in¬ 
sure that only the operating mode is common for both 
sections while the higher order mode spectra are different, 
i.e. the frequencies and Q-factors of the HOMs are differ¬ 
ent. Due to the resonant coupling between cells located 
on different axes, the voltage on each axis is nearly the 
same for the operating mode only, while it varies strongly 
for all non-overlapping higher order modes. Indeed, the 
cells act as a single structure at a common resonant fre- 
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quency (in our case it is the frequency of the operating 
mode) but they behave as separate structures at most 
other frequencies. This allows decoupling the accelerat¬ 
ing and decelerating structures for HOMs, thus breaking 
the positive feedback loop. 

Let us first define the model and the approach we use. 
When a bunch passes through a cavity it is decelerated 
by the self-induced voltage 


Vq = J^E{f,z,t)dl ( 1 ) 

where dl is defined along the beam trajectory and 
E{f,z,t) is the vector of the electric field seen by the 
beam. In the asymmetric system discussed the voltage 
induced is in general different for each axis, and will be re¬ 
ferred as Vq’^ with superscripts indicating the axis num¬ 
ber. The cavity geometry parameter R/Q which influ¬ 
ences its performance is presented as 


R/Qi,2 = 


ojeJylEir^l^dV 


( r \ 2m 

-) 

UJT J 


( 2 ) 


where m is the number of full wave azimuthal variations. 
The case where m — 0 corresponds to the monopole 
mode. The parameter R/Q varies from one axis to an¬ 
other as we are considering asymmetric cells. In a con¬ 
ventional symmetric system, which is either made of two 
identical sections, or the same cavity is used for acceler¬ 
ation and deceleration, the R/Q parameter is constant. 
Considering that the instantaneous cavity energy change 
is equal to the instantaneous bunch energy change, one 
can write (for the case of a single bunch inside the cavity): 


AUq = -AUc = 


/ c 1 

\ujr) 2uj R/Qj 


( 3 ) 


variations of shapes are small and leading only to the 
shift of eigenmodes frequency position while the uncou¬ 
pled modes’ transverse field structures will be considered 
to be the same. This also means that the same HOMs 
will be excited by the bunch inside of the accelerating 
and decelerating structures but at different frequencies 
and thus no positive feedback loop will be possible via 
the coupler (modes are fully separated and not overlap¬ 
ping). As the two structures have no feedback loop be¬ 
tween them, no beam instability (in conventional terms) 
will develop. However, those fields can still provide a 
large enough transverse kick to cause the bunch to hit 
the beam-pipe walls. In addition there is also a single¬ 
pass regenerative BBU where an instability develops in 
a single cavity. This takes place when a kick at the en¬ 
trance of the structure causes an offset at the exit of the 
structure, thus exciting a large dipole wake meaning that 
the following bunch will experience this wake field and its 
trajectory may deviate from the designed trajectory even 
further, leading to degradation of the energy recovery or 
even break-up of the bunch transportation. 
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where j = 1,2 indicates the number of the axis. For 
the particular case of the structure suggested, the cavity 
energy gain may vary if the modes’ eigen-fields at the 
axis 1 and 2 are different. The bunches moving along 
z will interact effectively with the modes having elec¬ 
tric field components collinear with the beam trajectory. 
The dipole modes with the transverse electric and mag¬ 
netic field components can contribute to the transverse 
momentum especially if there is transverse offset of the 
beam trajectory. The transverse momentum can also be 
gained via interacting with magnetic field of the HOM 
generated for instance by the previous bunch. If the 
bunch trajectory deviates from the designed trajectory 
which we will refer as r = \/= 0) the energy lost 
by the bunch to the HOMs will increase and the bunch 
will move further away from the axis. In this paper we 
are considering a small assymmetry of the cells, i.e. the 
cells on different axes have different shapes, however the 


FIG. 2. Schematic of asymmetric ERL structure with ac¬ 
celerated and decelerated bunches. The bunch without red 
filling indicated the bunch which deviated from the central 
line trajectory. 


In FIG. m the schematic of possible bunch trajectories 
are shown with r = rg = 0 representing the designed tra¬ 
jectory while Aa: deviation from the designed trajectory. 
A ‘transparent’ bunch in this figure illustrates a deviation 
of possible trajectory from the ‘design orbit’ (corresponds 
to a bunch motion along both central lines) which can be 
due to the excitation of high order modes or trajectory 
perturbation at the interaction point (IP). Further in the 
text, considering azimuthal symmetry of the system and 
to simplify notations (without losing generality) we con¬ 
sider only the xz plane (i.e. y = 0 and r = x) resulting 
in r{t + to) = x{t -\- to) = Rii^o -h Ri 20 . 
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A. Transverse deflection of a single bunch 

The deviation from the ‘design orbit’ can be calculated 
using 

a;(f + to) = Riixo + Ri20 (4) 


electrons passing through the accelerating gap L over 
time T. To evaluate Vj_ and thus the bunch deviation ac¬ 
curately, a full 3D eigenmode analysis is required. Here 
we will make few approximations to estimate a single 
bunch trajectory deviation. First of all, we take into ac¬ 
count monochromaticity of the field i.e. oc then (O 
can be rewritten as: 


where to is the time required for an electron to go from 
port 1 to port 2, 0 is the angle gained by electrons at 
the cavity exit, Rn and R 12 are the elements of the 
transport matrix that relates the angle at the cavity exit 
and bunch transverse position x at port 2. 

For simplicity but without losing generality we assume 
that at the entry point all bunches are having finite but 
negligible (in the first approximation) dimensions and 
move along the central axis. The finite transverse dimen¬ 
sion of the bunch will result in generation of transverse 
momentum leading to the shift of some electrons from 
the ‘design orbit’ of the system defined by the transport 
R matrix (we are ignoring for now some high order cor¬ 
rections associated with Rn term) 


a:(t-I-to) = R- 12 ^ (5) 

Here, for clarity only, we neglect the coupling of the 
offset and vertical/longitudinal beam parameters and as¬ 
sume that dl =dz. Taking into account that the beam is 
relativistic and its longitudinal velocity is close to speed 
of light, we can estimate the angle at which the bunch 
leaves the cavity as: 


Vj_ =- -e sin (j)T 

(jj 


f VxE^{r,(p,z)dz 

Jo 


(9) 


where = ojT /2 is a half transient phase. As we are 
looking at the maximum kick observed by the electrons 
it is clear that it can be achieved if 4>t = 7r/2, so the 
transverse potential becomes 


V± 




V±E^(r,(p,z)dz 


( 10 ) 


indicating that the effective transverse potential is 7r/2 
phase shifted from the longitudinal (accelerating) field. 
In the case of multipole fields, the following conven¬ 
tion is used to describe the longitudinal electric field: 
E^{r,ip,z) = l/2^{Ez{r, z)) cos{m(p). The polarisation 
axis is chosen such that the maximum electric field is at 
ip = 0. At this point a coordinate system change is also 
made to cartesian coordinates for simplicity and to match 
©. Making use of m and assuming (p = 0 such that the 
bunch see the maximum field excitation, the transverse 
potential can be written in terms of the longitudinal volt¬ 
age 


el4 ^ VI 

'^~W 


( 6 ) 


where (eVF) is the full bunch energy and V± is the effec¬ 
tive transverse voltage seen by the bunch. Using m and 
taking into account that xq = 0 one gets the expression 
for the radial deviation of the bunch from the second axis 


Ax — Ri 20 — Rl 2 -j^ (7) 

We note that if we use (y, z) plane it would transfer into 
Ay = R 340 — R34-^ 

The effective transverse potential which leads to the 
change of the total transverse momentum and thus tra¬ 
jectory deviation can be expressed as 

dt [ V±Ez{r,p,z,t)dz (8) 

Jo 

where = ^r+ ^^cp. 

Here for convience we keep cylindrical coordinates 
however if inserted into ([7]) a coordinate change may be 
required. This introduced potential is observed by the 



Vx 



( 11 ) 


We also assumed that a bunch propagating along the 
‘design orbit’ is perfectly timed to ‘see’ both the maxi¬ 
mum accelerating and decelerating potentials. Any tra¬ 
jectory perturbation due to interaction at IP is ignored. 
Assuming that bunches arrive to port I on axis, the tra¬ 
jectory shift at the second port is given by: 




Ax = Ri 2 - 


fo ^E,{x,z)dz 


W 


( 12 ) 


The expression (fTTl) gives an opportunity to make a 
first estimation of an upper bound for the maximum R 12 
parameter knowing some basic properties of the cavity 
and assuming that E^{x, z) oc cos{kzz) (pill-box like cav¬ 
ity). Evaluating the integral in ((T7|) we get: 



w 





( 13 ) 


where = kzL is the electron bunch phase shift at the 
cavity exit and kz = \/{ojjcY — k^. We assume that 
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fc;; = 0 as in this case (fiz = kzL = 0 and function sincc/x 
in m has maximum at this value of x. Thus, for all 
other values of kz, the function will be smaller and Ax 
will takes values from 0 to Ax as defined in m 


2c d L 

Ax^- ■7^V{x)—B.i2 (14 

w ax IV 

This expression shows the link between bunch’s trajec¬ 
tory deviation and machine’s parameters. Therefore, to 
reduce the deviation one may either reduce the parame- 
ter or limit the R 12 parameter for a given aperture 
(diameter D 2 ) of deceleration cavity (port 2 ) as shown 


Ri2 1 ) 2 ^ 


2fc.L2 |j.y| 


(15) 


I dx 


One can further evaluate the expression by assum¬ 
ing polynomial representation 0 of the potential, i.e. 

M 

V = where m = 0 is associated with the 

m—0 

monopole mode and m = 1 with the dipole mode (we 
ignore the other modes based on the assumption of small 
deviation from the ‘design trajectory’). Let us recall that 
we assumed r = x. Substituting this into m one gets 


Ri2 ^ D 2 W 


1 




M 

E 

m—1 


^m+1 


(16) 


while the time deviation from an optimum time along the 
‘design trajectory’ is 


^ 1 /Ri20Vti Ax Ax 


(18) 


Multiplying m by the frequency of the operating 
mode ojQ and taking into account that oJoTg = 27rn where 
n is the number of field oscillations in the cavity during 
the bunch travel from the accelerating to the decelerating 
arm, (US can be transformed: 


A^.= (^)%n (19) 

It is clear that if there is no deviation there is no phase 
variation from the optimal value. Noting that n has the 
meaning of the number of wavelength of operating modes 
Ao = 27rc/a;o along the optimal path Sq, we can present 
(IT^ in a slightly different form 


A(j) 


(Ri20)Ax 

SoXo 


( 20 ) 


Taking into account that Ac/) <C tt, one can write a 
second condition for R 12 in order to guarantee a small 
deviation of Ac/) and thus effective deceleration. We also 
note that to avoid beam dumping before deceleration cav¬ 
ity the condition Ax < D needs to be satisfied and thus 
we can rewrite (EOl) as 


It is clear that this condition is not sufficient for the 
ERL to operate as the bunch with a trajectory deviated 
from the ‘design orbit’ will be delayed due to the longer 
path given by the transport matrix element R 52 and its 
deceleration (energy recovery) will be affected, leading 
to possible interruption of ERL operation even if the 
bunch passed through the decelerating section. If the 
beam enters the decelerating cavity with a deviation of 
the transverse position xq + Ax then the transient time 
Tg from acceleration to deceleration sections is changed 
as (Tg-I-Ar) and deceleration will shift from its optimum. 
This phenomenon is a bunch de-phasing and knowing Ax 
one may estimate the bunch time deviation At from the 
time travel along the design trajectory Tg due to trajec¬ 
tory deviation Ax. Ignoring transport matrix elements 
such as R 52 but keeping R 12 , i.e. neglecting any beta¬ 
tron oscillations in the transport line and assuming that 
6x So, where is the length of the design trajec¬ 
tory, and that the electron bunch travels with the speed 
of light c, one finds that the total length of the deviated 
trajectory is 


5 = x /52 + (Ri 20)2 ^ 5o 


1-f 


1 f 


(17) 


A(/) 


(Ri 26 ')Ax 

^oAo 




TT <C TT 


( 21 ) 


In conventional accelerators condition (HU is 
“stronger” and thus sufficient. The effective trans¬ 
verse potential due to multiple bunches can also be 
derived. The derivation can follow the approach 
presented in [l^ and can be found in the Appendix El 


III. HOMS OVERLAP AND BBU START 
CURRENT 

In this section we will employ RLC circuit approach to 
describe the system. Let us now consider that in reality 
the ROMs of the accelerating and decelerating sections 
will overlap due to their finite bandwidths. In FIG. [3] a 
schematic of such an overlap is shown. At the shaded 
(overlap section) part of the spectra one may expect a 
leakage of the energy from one arm of the cavity into 
another. This will lead to a coupling of energy between 
the accelerating and decelerating sections, providing a 
feedback path for instabilities to develop. One such in¬ 
stability is multi-pass regenerative BBU. 

Here we meet our first issue when trying to describe 
a dual axis machine. The standard RLC equations are 
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FIG. 3. An example of two modes with different frequencies 
which overlap due to their finite bandwidth. The shaded re¬ 
gion is helping illustrate where two modes overlap but does 
not indicate the magnitude of coupling. 


not valid as the impedance is dependent on which cavity 
axis the beam traverses. However, for a given energy and 
for a given mode the ratio between the voltages on each 
axis is constant, a scaling factor can be introduced. In 
essence the system acts like a transformer with a load 
and a current source at both the primary and secondary 
axis, with a ratio of Nsp. Considering the voltage at the 
accelerating section Vi and at the decelerating section V 2 
for a given stored energy we get 


C 2 = NspVi (22) 


ing into account that the shunt impedance of the uncou¬ 
pled acceleration and deceleration sections are the same, 
Ri = R 2 , this leads to the shunt impedance of each struc¬ 
ture, if coupled, to be different when seen from the other 
axis. Indeed, the impedance of the secondary structure 
(located on the 2nd axis) when viewed from the source 
on the hrst axis is i? 2 , and can be written as [Tsil 


R '2 = RiNIp = R 2 /NIP ( 24 ) 

Also the beam across one accelerating/decelerating gap 
can be looked at as being a current source in parallel with 
both impedances i?i _2 (loads). In this case the beam will 
apply voltages across the gaps while the stored energy, 
and the losses will be also split between the both loads. 
Such a description of the system allows one to define 
the total impedance of the whole system as a sum of 
two parallel impedances and find the voltage on the first 
axis, Vi, and second axis, V 2 , induced due to the bunch 
propagation along the first/second axes. The expression 
for Vi can be calculated from (1^ and is given in this 
case as: 


Vi=h 


Ri 

1 + NIp 


( 25 ) 


And the voltage across the decelerating gap due to the 
beam in the accelerating gap is similarly given by 


Cavity 1 Coupling cell Cavity 2 



V2 = NspV' = NspVi = W2 ( 26 ) 

1 - 1 - Ngp 

The total effective transverse impedances for both gaps 
(combined in parallel) R±i and i?_L 2 can hence be intro¬ 
duced for the accelerating and decelerating gaps respec¬ 
tively 


R±i = 


Ri 

1 + NIp 


( 27 ) 


FIG. 4. A circuit diagram of the dual axis structure. 

From the primary current source the voltage of the 
secondary is stepped either up or down depending on the 
transformer ratio. In this context, it is useful to define 
the voltage in the decelerating section as seen from the 
accelerating section such that we can view the system as 
two parallel circuits. The voltage on the second axis, V 2 , 
when looked at from the primary axis (i.e. the bunch 
propagates along first axis) is transformed to which 
will be the same as Vi and therefore 


Vi = V/ = 


Nsp 


( 23 ) 


First lets calculate the shunt impedance of each struc¬ 
ture, ignoring the losses in the other structure. Tak- 


R±2 = NspRp, = ^^j;^ ( 28 ) 

Let us apply these expressions to the calculation of 
multi-pass regenerative BBU. A bunch traversing the ac¬ 
celerating gap of the cavity will experience a transverse 
force due to the previous bunches excitation of a dipole 
HOM and will deposit energy into the cavity. It will expe¬ 
rience an acceleration or deceleration dependent on the 
offset of an individual bunch from the electrical centre 
of the cavity, Xa- The sum wakefield due to the previous 
bunches (we are looking at a single dipole HOM) has am¬ 
plitude Vd, phase (j) and the accumulated energy change 
of the bunch in terms of voltage is AC/i, 


AC/i = qXaVd cos 4> 


( 29 ) 
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Let us note that the phase (j) is between the bunch cen¬ 
troid and the peak accelerating EM field excited by the 
previous bunches The trailing electron bunch will have 
a deflection which is defined by an effective transverse 
voltage (as discussed above) and is given in this case as: 

V^i =-Vdsu\(j) (30) 

UJ 

We note that the transverse effective potential and de¬ 
celerating voltage are —7r/2 shifted in respect to each 
other (term -i in front of the expression) as given by 
Panofsky-Wenzel theorem 0- A bunch deflected by this 
effective potential will arrive at the decelerating gap with 
a transverse offset from the design orbit, Ax, which (as 
we discussed) can be calculated for a given beam energy 
W, using the transfer matrix R, 


Ax = Ri 2—sin^ (31) 

Similarly as we done above, we assumed that RnXa 
is small and therefore, neglect this term. Since the two 
gaps are coupled, the dipole mode also has a longitudi¬ 
nal voltage component in the decelerating gap leading to 
acceleration or deceleration of the bunch, and the max¬ 
imum energy A172) deposited in the mode by the beam 
traversing the decelerating gap, is given by 


A172 = q{xd + Ax)NspVdC{Tg) (32) 

where C[Tg) = cos{(j) + i^Tg), Xd is the distance between 
the HOM electrical centre and the design orbit in the de¬ 
celerating axis (in general Xa ^ Xd)-, Tg is the time travel 
from the accelerating to the decelerating cavity. We as¬ 
sume that the decay of the field in the cavity between 
the bunches going through each axis is small. For a pure 
dipole mode, if Vd grows in time then so does the bunch 
offset Ax, and hence to calculate the BBU threshold cur¬ 
rent we only need to work out the power balance in order 
to derive the condition for the instabilities not to grow. 
The average power over one RF cycle (assuming every 
bucket is filled) is 


hVdXa cos (j) + Io{xd + Ax)NspVdC{Tg) - Pc <0 (33) 

where Iq is the beam current, Pc is the energy losses 
due to the finite Q-factor of the cavity. The first term 
on the left hand-side of the expression indicates the EM 
field energy gain inside the accelerating section, while 
the second term shows the energy change of the EM field 
inside the deceleration section. It should be noted that 
for the case of every bucket being filled the 1st and 2nd 
terms may not necessarily be due to the same bunch. 
Inserting the bunch offset at the second axis due to the 
kick on the first axis we obtain 


loVdXo cos (j)+ 


loNspVd 


Xd 


Ri2tt 7—sm0 
W UJ 


C{Tg) - Pe > 0 


(34) 


The cavity power losses for a given load are defined by 
Ql (Q-factor which takes into account both internal and 
external losses due to HOM couplers) which is given by 


where R/Qpi is the transverse R/Q for a dipole mode 
(as defined above). Inserting this into equation (IMl) and 
substituting for Vj_ yields 


loVdXa cos cj) 

+ loNspyd 


vl 


Xd + R-12 


> 0 


Vd 

w 


0 sin<)> 




(36) 


RlQxiQh 

Rearranging for Jq yields 


do > 


cVj_ 

ujRjQpiQh (xa cos c/) + NspxC(Tg)) 


where 


(37) 


Xd , Ri2 . , 

We note that, since Vj_ and (f> are not known, in order to 
find the solution, a full wakefield calculation and analysis 
are required in such an asymmetric structure. However, 
if we consider (just for our case) that both of the offsets 
on the first axis (i.e. bunch centroid offset and the offset 
between the electrical centre and the bunch design orbit) 
are small as compared to the bunch offset at the second 
axis, then due to the transverse kick near the BBU limit 
the non-equilibrium current can be simplified to 


cW 

NspujRi2R/Q±iQLC{Tg) 


(39) 


This simplification does however mean that the start 
current calculated is an over estimate by a factor of ap¬ 
proximately {xd + a)/a where a is the beam aperture. 
This equation is however still a function of the initial 
transition phases (j) and ujTg. Let us look at the worst 
case (i.e. the largest transverse kick at the first axis when 
(j) = irj 2 ), which leads to the following estimate 


cW 

do > 


AspujRi2R/Q±iQl cos(ujTg) 


( 40 ) 












If ojTg = 27rm and cj) = n/4 the equation (l40l) becomes 


NspujR/Q±iQLii^i2 

Expanding this expression using the transformer ratio 
and the relation between the impedances of a single sep¬ 
arate structure (as discussed above and not the effective 
transverse impedance of the whole cavity) one will get 


done to verify the results observed and illustrate some 
of the conclusions. For all cases, curvilinear tetrahedral 
mesh-elements were used. We have also investigated the 
Wakefield generation. The wakefield was Fourier analysed 
and compared with location of eigenmodes of the ‘cold’ 
system. Taking into account the repetition rates of the 
system and considering analytics developed, we show the 
possibility to increase further the BBU start current via 
the introduction of the additional losses for HOMs. 


(1 + NIp) cW 
NspUj'R.j^2R/Q±iQl 


(42) 


Comparing it with the start current of symmetric cav¬ 
ity, Nsp = 1, one notes that the start current of an 
asymmetric cavity is increased. In a symmetric cavity 
each HOM has equal field in each structure {Nsp = I) 
leading to the following relation between the currents: 


I, 


asymmetric 


> 


(1 


^S p) ^symmetric 


2Nsp 


(43) 


Analysing the expression, one notes, that the ratio be¬ 
tween the currents has the minimum value of 1 if Nsp = 1 
(a symmetric cavity). This can be interpreted in the fol¬ 
lowing way. If Nsp < I, the energy feedback from the 
second cavity is smaller, leading to start current increase 
(as the denominator is proportional to Nsp). On the 
other hand, if Nsp is large, the kick on the first axis be¬ 
comes small and to get BBU, the start current should be 
increased, which manifests itself in the Ngp term. One 
also notes that for Nsp = 1 the cavity has twice the 
power losses of a standard cavity and thus the start cur¬ 
rent is doubled compared to a recirculating machine. 


A. Cavity Design 

Let us consider the design of the ERL structure pro¬ 
posed. The structure is created from two typical ellip¬ 
tical SRF cavities with different parameters joined by a 
coupling cell. A typical elliptical SRF cavity is formed 
from a number N of mid-cells plus end-cells to take into 
account the coupling to the beam pipe. For the asym¬ 
metric structure, the design is more complex as shown in 
FIG. 0 A key step in designing the structure is to make 
the mid-cells for each axis different but still maintaining 
the same operating mode. 


Coupling cell Mid cells 



FIG. 5. Schematic of one axis of the structure. 


IV. NUMERICAL MODEL OF WAKE FIELD 
GENERATION IN ASYMMETRIC CAVITY 

A bunch propagating inside the structure excites wake 
fields which can be described as superposition of eigen¬ 
modes of the cavity i.e. E{r,z) = where Cg 

represents a coefficient and Eg corresponds to the elec¬ 
tric field for the eigenmode with indice s. The excitation 
of a single eigenmode is discussed above and calculation 
of complex wake field considering all modes is outside 
the scope of analytical theory and can be done using nu¬ 
merical modeling. The goal of the numerical modeling 
in this case is to identify the most dangerous modes in 
the system described as close to the real case as possible, 
using eigenmode solvers, and using analytics developed 
to suggest the methods to suppress these modes. The 
simulations were performed using the ACE3P electro¬ 
magnetic suite developed at SLAG. Taking into account 
that such systems were not modelled previously we con¬ 
ducted investigations using both frequency (Omega3P 
[13) and time (T3P [l^) domain approaches. This was 


The design process involved using the mid cell design 
from a conventional cavity (in this case, the TESLA cav¬ 
ity shape, operating frequency of 1.3 GHz) and to 
vary the parameters, finding two shapes which had the 
biggest difference in geometry yet still shared the same 
operating mode. The mid-cells are required to have the 
same frequency and R/Q to ensure equal voltages in both 
modes. The R/Q does not have to be exactly the same 
but the frequency must be. A change in R/Q will re¬ 
sult in different field amplitudes in the two sections. The 
length of the half-cell is kept to a quarter of a wave¬ 
length such that the designs can be scaled to other fre¬ 
quencies. For example, if the desired frequency was 1.95 
GHz, each parameter can be multiplied by a factor of 
1.3/1.95. Additional tuning is needed for the end-cells 
where the length is no longer half a wavelength. Of im¬ 
portance is to try to separate the higher order mode spec¬ 
trum as much as possible. Once two shapes were found, 
a dispersion diagram was created for each design and an 
example is shown in FIG. [71 Each line on the disper¬ 
sion represents the passband for a particular mode in the 
case of an infinitely long periodic structure. A passband 
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shows the frequency range in which a mode can propa¬ 
gate with a particular phase advance between adjacent 
cells in a structure. When the structure is used with a 
finite number of cells, a pass band will split into a finite 
number of modes that lie on the line. We assume that 
the optimum scenario is that the pass bands will have a 
frequency separation on the order of one or a few MHz. 


Midcell 


Coupling cell 



FIG. 6. Schematic of the mid cell (left) and the coupling cell 
(right). 


2.4r 


,_^ 2.2 

K 2.0 
O . 
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a> 
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0) 


Design 1 
Design 2 


7r/5 


27r/5 37r/5 

Phase 


47r/5 


FIG. 7. Dispersion diagram for the two mid cell designs. A 
clear separation between the passbands can be seen with the 
exception of the first passband in which the operating mode 
is shared. For a symmetric structure the passbands would be 
identical. The dotted line shows the speed of light curve. 


Parameter Axis 1 cell [mm] 

Axis 2 cell [mm] 

Mid cells 

Req 

103.3 

103.3 

A 

42 

42 

B 

42 

43.1 

Riris 

35.75 

37 

a 

12.75 

11.75 

b 

18 

20 

1 

57.7 

57.7 

End cells 

Req 

103.3 

104.3 

A 

42 

42 

B 

42 

43 

Riris 

39 

39 

a 

12.75 

11.75 

b 

18 

20 

1 

58.54 

60.96 

Mid to coupling cells 

Req 

103.3 

104.3 

A 

42 

42 

B 

43.4 

43.5 

Riris 

35 

35 

a 

12.75 

9.69 

b 

18 

20 

1 

57.7 

57.7 

Goupling cells 

A 

48.052 

48.052 

B 

29 

29 

Riris 

35 

35 

a 

9.6 

9.6 

b 

10.152 

10.152 

l 

57.652 

57.652 

Ls 

150 

150 

Lb 

111 

111 

End coupling cells 

A 

47.5 

47.5 

B 

29.76 

29.76 

Riris 

39 

39 

a 

9.945 

9.945 

b 

9.945 

9.945 

1 

57.652 

57.652 

Ls 

150 

150 

Lb 

111 

111 


TABLE I. Parameters used to construct numerical model. 


The design of the coupling cell is more complex. Its 
shape is racetrack-like in order to couple both accelerat¬ 
ing and decelerating sections with the axes of the sections 
separated by distance 2Lf, (FIG. [b]). The wall angle be¬ 
tween the two ellipses is forced to be 90°. The example 
is shown in FIG. [5] and TABLE H] The main purpose of 
the design is to couple strongly the sections at operat¬ 
ing frequency and effectively split them at all the other 
frequencies at which HOMs are located. 

To create the end cell, the parameters of the mid cell 
were modified to take into account the beam pipe di¬ 
mensions. Another change in dimensions is also needed 
to join the mid cell to the coupling cell. The mid cells for 
each axis have a different iris radius, however the cou¬ 
pling cell is chosen to be symmetric for simplicity and 


thus another cell is needed to match the iris radius with 
the coupling cell. The example of such dimensions is 
shown in TABLE HI One notes that these dimensions can 
vary depending on the operating mode frequency. 

B. Numerical analysis - Frequency Domain 

As we mentioned in the introduction to this chapter 
both types of analysis i.e. frequency and time domain 
were performed (for verification and illustration purpose) 
to study these structures. FIG. [8] shows the desired op¬ 
erating mode (at 1.3 GHz) transverse structure together 
with the Ez field flatness along the longitudinal axis z 
(FIG. |S]). The field flatness was calculated ignoring the 
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coupling cell and for both axes is better than 98% which 
is acceptable in most practical cases and can be further 
improved if required. It should be noted that the design 
for this paper provides a proof of principle and is by no 
means a fully optimised structure. 

High efficiency SRF ERLs require very high quality 
factor of the fundamental mode. Thus, there will be a 
very strong sensitivity to vibration and pressure known as 
’’micro-phonics”. To the operating mode the two cavities 
look like a single eleven-cell cavities rather than two cou¬ 
pled cavities, as the coupling cell is also resonant. This 
is no different than a normal nine cell cavity where the 
microphonics can shift the cell frequencies independently. 
Tuning becomes more difficult as the cavity gets longer 
as the modal frequencies in the passband come closer to¬ 
gether, which is why the design shown here is limited to 
five cells on each side and strong coupling is used. 

Due to the different cell designs the amplitude of the 
electric field in each axis is slightly different (FIG. 6b) 
which can also be beneficial (i.e. possibility to vary the 
field in accelerating and decelerating cells) for energy re¬ 
covery after beam interaction at the IP. We have calcu¬ 
lated and analysed the first 100 eigen-modes of the sys¬ 
tem. For each mode, the complex longitudinal voltage 
I/|| jj(r) was extracted in order to calculate the longitudi¬ 
nal R/Q and transverse R/Q for each axis.The longitu¬ 
dinal voltage was extracted using 



while the longitudinal R/Q was calculated using 


R/Qn 


^nUn 


(45) 


Considering the fact that the voltage for a dipole mode 
scales linearly with x, the transverse voltage Vj_,n can be 
written using the longitudinal voltage as 


= [f)|.n(r)-k)|.„(0)] (46) 

UJn'^ 

In a similar way to the longitudinal R/Q, the trans¬ 
verse R/Qx,n can be defined as 


R/Qx,u = (47) 

The transverse R/Q was evaluated at 1 mm offsets in 
the horizontal (x) and vertical (y) directions from the 
centre of each axis. FIG.Oshows that the i?/Q’s are not 
equal for each axis as desired. A cluster of modes is visi¬ 
ble between 1 and 1.4 GHz. The majority of these modes 
lie inside the fundamental pass-band (around 1.3 GHz) 
and include the operating mode. On either side of this 


pass-band two modes have equal R/Q (at ^ 1.1 GHz 
and ^ 1.5 GHz) are present. However these modes are 
conhned/localised to the coupling cell and damping is 
thought to be possible. This will be illustrated and dis¬ 
cussed below. The next cluster of modes has high trans¬ 
verse R/Qx suggesting these are dipole modes. These 
dipole modes having the highest R/Q± are of most con¬ 
cern for preventing the development of BBU. 



FIG. 8. A contour plot of the electric field distribution for 
the operating mode (left) and the electric field along each 
axis (right). 
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FIG. 9. The longitudinal and transverse R/Q evaluated along 
both axes. 

Looking at the worst scenario we note that the trans¬ 
verse R/Q oi one of the mode has an R/Q of 72 H. In 
conventional system this mode would have the same R/Q 
along the second axis leading to BBU development. How¬ 
ever in this case the R/Q along the second axis is only 
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Axis 1 Axis 2 Axis 1 

Axis 2 Axis 1 

Axis 2 

Frequency 

R/Q 

R/Q±,x 

R/Q±,y 

GHz 

n 

n 

n 

Highest R/Q 

1.3 

348.71 

301.51 

0.0675 

0.0365 

0.0074 

0.0 

1.29943 

231.71 

247.59 

0.0014 

0.0059 

0.0003 

0.0048 

1.09966 

32.622 

32.367 

9.5769 

9.0660 

0.0166 

0.0055 

1.29532 

21.075 

23.878 

0.0014 

0.0267 

0.0281 

0.0333 

1.48554 

20.337 

20.429 

12.094 

12.360 

0.0026 

0.0001 

Highest R/Q±,x 

1.70216 

0.0035 

0.0127 

65.207 

0.8680 

0.0134 

0.0004 

1.74343 

0.0211 

0.0069 

61.997 

0.4792 

0.0294 

3.8679 

1.87193 

0.0050 

0.0029 

35.500 

0.0810 

0.0555 

0.0002 

1.85436 

0.0181 

0.0091 

17.329 

0.3260 

0.0208 

4.2119 

1.48554 

20.337 

20.429 

12.094 

12.360 

0.0026 

0.0001 



Highest R/Q 




1.73192 

5.4419 

1.4736 

0.0005 

0.0001 

72.089 

0.3764 

1.68526 

1.6890 

9.9178 

0.0024 

1.8274 

36.312 

0.6537 

1.78142 

1.5499 

8.7076 

0.0039 

0.0070 

25.636 

0.1329 

1.87103 

1.6368 

7.9211 

4.8525 

0.0037 

22.491 

4.0005 

1.8523 

7.7902 

6.4131 

0.0033 

0.0001 

15.388 

0.1740 


Accelerating/decelerating cells 

I I _ I I Coupling cell 

mV 

Beam pipe 


FIG. 11. A mode confined to the coupling cell which is sym¬ 
metric. 


alternative is to perform a time domain calculation in 
which the beam pipe boundaries allow waves above the 
cutoff frequency of the beam pipe to propagate. The 
fields generated by a charge on a trailing particle can 
also be easily calculated and compared/verified with the 
data observed from frequency domain analysis. 

To study wake-field excitation let us consider a Gaus¬ 
sian bunch with charge density 

1 (s-3o)^ 

pis) = (48) 

VZttct 


TABLE II. The modes with the top five highest RfQ, RIQ±,a: 
and R/Q±,y. 

0.38 n. This confirms that due to the asymmetric design 
of the system it works as expected, namely confining the 
HOM to just one of the axis. FIG. [TUI shows a transverse 
and longitudinal slice of this mode. 



FIG. 10. A transverse and longitudinal slice showing the elec¬ 
tric field contour plots for a mode at 1.73 GHz. The transverse 
slice also shows the magnetic field as indicated by the cones. 

As mentioned above and shown in FIG. [U] there are 
modes (due to symmetric shape of the coupling cell) 
which have the same R/Q along each axis. These modes 
are confined to this cell (FIG. [IT]) resulting in the R/Q 
being equal along each axis but also meaning that such 
modes can be easily dealt with by applying an appropri¬ 
ate coupler or an absorber. An example of such a mode 
(located at frequency 1.4855 GHz) is shown in FIG. [TT] 
and it is clear that locating a coupler at the middle of 
the coupling cell (for instance) will not affect the beam 
transportation and will effectively damp this mode. 

C. Numerical analysis - Time Domain 

One of the limiting factors of frequency domain cal¬ 
culations is the use of closed boundary conditions. An 


where a is the bunch length (we assumed that cr = 1 cm), 
s is the distance along the path. The bunch is then 
launched along one axis of the system with a vertical 
offset assumed to be 1 mm. By this way the full 3D 
model is simulated meaning both monopole and dipole 
like modes are excited. The simulations were repeated 
with the bunch traveling along the second axis in the op¬ 
posite direction. The longitudinal wakefield was calcu¬ 
lated using Weiland’s indirect scheme and the wake- 
field potential was Fourier transformed/analysed to ob¬ 
tain the complex impedance spectrum using 

I 

Z||(a;) = —-/ IF||e-“^/^ds (49) 

cp(a;) 

The Panofski-Wenzel theorem [l^ was then applied in 
order to calculate the transverse impedance spectrum 

-Zi(cc) = Vj.^||(w) (50) 

c " 

where V_l = + -^y. The schematics of simulations 

discussed is shown in FIG. m The two simulations are 
similar but in the first simulation, bunch 1 is launched 
along axis 1 in the positive Z direction and the wake is 
calculated along axis 1 and 2 with a frequency resolution 
of 2 MHz. The simulation is then repeated (model 2) 
for the opposite case (with the same resolution) where 
bunch 2 is launched along axis 2 in the negative Z di¬ 
rection. In both models a time-step of 2 ps was used 
to observe the convergence of the results. The spec¬ 
tra of longitudinal Zyy and transverse Z±ij impedance 
observed from these simulations are shown in FIG. m 
The subscripts ij indicate in both cases the number of 
axis and the bunch as discussed. Let us note a good 
agreement between FIG lUI and (TUI indicating a good con¬ 
vergence of the results observed from the time and the 
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frequency domain approaches. One can easily identify 
different families of modes and see a good correlation be¬ 
tween the amplitude of impedances and R/Q parameters 
of the modes. It is also interesting to note that the spec¬ 
tra of the impedances of the wakefields associated with 
HOMs (frequencies above 1.5 GHz) measured on differ¬ 
ent axes are different (light and dark blue lines can be 
easily distinguished in the figure). 



FIG. 12. A schematic of the simulations setup. 




Indeed it can be seen from FIG. M that the operat¬ 
ing mode’s impedances (1.3 GHz) are roughly equal. A 
number of other modes can been seen with almost equal 
coupling to the beam including a coupling cell modes 
(equidistantly located on either side of the operating 
mode) and a band starting at 2 GHz. The last group 
of the modes is also associated with the symmetric cou¬ 
pling cell and can be easily damped if necessary. A band 
of cavity modes which are of our concern is located be¬ 
tween 1.6 and 1.9 GHz. Their impedances spectra vary as 
desired and thus they may not be considered as a threat 
to the system stability. Also these modes are dipole like 
modes and thus their contribution is better shown in the 
transverse impedance spectrum in FIG 1131 

The transverse impedance spectrum still shows the op¬ 
erating mode but to a much smaller degree. The effect of 
the cavity modes between 1.6-1.9 GHz is now much more 
visible. As desired, the frequency and coupling of these 
modes are different depending on the axis. There is now 
a lower order mode at 1.1 GHz which is the monopole 
mode of the coupling cell. It has a large impedance, is 
symmetric and as the beampipes are not at the centre, it 
will have a transverse component. 


V. CONCLUSIONS 


FIG. 13. The longitudinal impedance (top) and transverse 
impedance (bottom) spectrum. 


The BBU start current can be calculated for each mode 
in the cavity, and the lowest value is taken as the start 
current for the machine. Let us first calculate the worst 
case {(j) = 'K jT) and uiTg = mr -\- it 12 where n is an in¬ 
teger. In FIG [14] the comparison of the start currents 
for symmetric and asymmetric cavities is shown using 
ujTg = 27rn. As we are only concerned with the modes 
with the lowest start currents we only show modes with 
start currents below 6 Amps. As can be seen the start 
current increases by a factor of 4. 

6 


5 


._,4 

t: 

^ 2 

I 

0 
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# Symmetric 
X Asymmetric 


1.4 1.6 1.8 2.0 

Frequency [GHz] 


In order to evaluate the benefits of such an asymmetric 
cavity we calculate the start current for a sample machine 
with an injection energy of 5 MeV and a cavity voltage 
of 5 MeV. We will assume an R 12 = 1 m. It is neces¬ 
sary to calculate both the BBU start current (calculated 
in section lYI) as well as current at which the wakefield 
drives the bunch into a wall due to the electrical centre 
shift of the dipole mode (as calculated in section uni. 


FIG. 14. Gomparison of the start currents for symmetric and 
asymmetric cavities is shown using Tg = 2?™ 

As a next step, one can look at the case where Tg is 
fixed. In order to maximize the energy recovery process, 
the delay time must be an integer multiple of the RF 
period at the operating frequency. Here we use the 10th 
sub-harmonic of 1.3 GHz giving a delay of 7.69 ns which is 
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FIG. 15. Comparison of the start currents for symmetric and 
asymmetric cavities is shown using Tg = 7.69 ns. 
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a path of just over 2 m. In FIG[Tn]the comparison of the 
start currents for symmetric and asymmetric cavities is 
shown for such a delay Tg=7.69 ns. The start current ob¬ 
served for the asymmetric machine of 3.8 A is dominated 
by the symmetric dipole mode at 1.72 GHz (eigenmode 
of the symmetric coupling cell discussed in section 4), 
limiting the benefits of the asymmetric cavity. However, 
this mode has a smaller R/Q than the operating dipole 
modes and also this mode can be effectively dealt with 
by either providing strong HOM damping or by alter¬ 
ing the cells such that the sum wake has a phase which 
is an integer multiple of tt. These methods are viable 
for the asymmetric system as only a single mode needs 
to be considered. As can be seen, the start current for 
all other modes is significantly higher for an asymmet¬ 
ric machine (3.8 A) compared to a symmetric machine 
(593 mA) by a factor of 5. One may argue that due to 
the asymmetry of the cavity the dipole modes electri¬ 
cal centre doesn’t align with the design orbit leading to 
the possibility of the induced wakefield to become large 
enough to drive the bunch into a wall even without a 
feedback path. However, using the equations derived in 
this paper one finds that the current at which the wake 
kicks the beam to a significant offset to dump the beam. 


i.e. 10 cm, is over 1000 A, and hence is not a limit to the 
machines maximum operating current. 

If the coupling cell were to be redesigned it could be 
possible for the start current to tend to infinity. One 
could design a coupling cell which would make use of 
magnetic coupling. However, a magnetic coupling cell 
would require a weld at a high current region and such 
approach was discounted, at present, for this design. For 
the design presented here, the coupling cell is also sym¬ 
metric. An antisymmetric coupling cell could also be 
designed to increase the start current threshold further. 

If the transit time is ignored, the limiting factor is 
driven by symmetric coupler modes however, in the case 
when the transit time is included, the length of the arc 
connecting the acceleration and deceleration arm can be 
chosen such the beam traverses the deceleration arm in 
anti-phase with respect to a symmetric mode, thus avoid¬ 
ing BBU. The limiting factor then arrises from the fact 
the Nsp is non-zero and thus, the BBU limit will be 
raised compared to a symmetric structure but not infi¬ 
nite. The reason Nsp is non-zero is because both cavities 
sections are coupled, eigenmodes will be of the full sys¬ 
tem and not confined to one arm. Modes that have high 
Nsp are probably caused by two modes very close in 
frequency space such that there is some coupling. 

To sum up: in this paper we introduced a new con¬ 
cept of asymmetric system of accelerating and decelerat¬ 
ing cavities coupled by a resonant coupler for ERL. We 
considered both numerical and analytical approaches and 
demonstrated the operational principles of new system. 
We have compared it with conventional symmetric de¬ 
sign and discussed the advantages of the asymmetric lay¬ 
out. Indeed it was demonstrated that asymmetry allows 
transporting through the ERL an electron bunch having 
significantly increased (at least by factor of 5) current 
without limiting BBU. This potentially is very attractive 
for applications in which it is important to increase THz 
and x-ray photon yield as well as bunch energy recovery. 
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Appendix A: Effective transverse potential due to 
multiple bunch interaction 

From the previous discussion it is clear that bunch de- 
hection and time delay are affected by the effective trans¬ 
verse potential which we discussed for a single bunch in 
the sections above. Therefore, to calculate Ax of bunch 
N-l-1 after N bunches passed through the system or the 
bunch delay at the decelerating section, the transverse 
kick due to held accumulation inside the cells needs to 
be estimated. We will note that are the func¬ 

tions which dehne TMon,— like and TEo„— like modes’ 
transverse structures. To derive the transverse dehec- 
tion, expression m will be used and we will consider 
the worst scenario i.e. maximum transverse kick. We will 
also consider HOMs of a single un-coupled cavity with the 
strongest contribution to the 6. This is arisen from the 
fact that the asymmetric structure is made of two cav¬ 
ities which are weakly coupled, away from the resonant 
frequency of the coupler. For the purpose of discussion, 
we use only the hrst few dipoles HOMs which have the 
highest impedance and hence dominates the long range 
wake allowing us to ignore the rest of the high frequency 
HOMs. We note that in the model described, after energy 
recuperation stage, the bunch is dumped into a collector. 
Due to this the system has the following specific times: 
repetition rate Trep , time decay T^ec of HOM (we as¬ 
sume that Q-factors of respective HOMs are the same 
in both cavities), and the bunch transit time, to between 
the accelerating and decelerating arms. Let us now define 
Tdec > Trep and Tdec > to while taking into account that 
all these times are much larger as compared with bunch 
transient time T through the structure. We note that a 
field amplitude seen by the following (second) bunch will 
decrease in time (dissipation losses due to finite value of 
quality factor Q for HOMs): 


l^2| = 




(Al) 


Vi ,2 are the voltages generated (3) by the first (driver) 
and seen by the witness bunch respectively and If Tdec = 


(A2) 

Here we use cylindrical coordinates to describe fields only, 
Ez is generated by the bunches and seen by the witness 
bunch. In (IA2I) the introduction of ‘discrete’ time for 
EM field accumulation i.e. proportional to Trep is possi¬ 
ble due to very short transient time through the struc¬ 
ture Trep T. To analyse the field accumulation we 
can introduce a decay a and detuning Sujn parameters 
a = l/Tdec, to = 2TTh/Trep + Sujn where h is an integer 
number indicating a number of full oscillation which EM 
field does at the operating frequency before next bunch 
enters the cavity. Clearly h is not connected to n (num¬ 
ber of bunches) and can be either arbitrary large or small 
(large or h —>■ oo is the case of single bunch approxima¬ 
tion) . Taking only a real part of the sum 

N 

sion can be rewritten as: X) ® 

n—1 

leading to the expression 


a~n'Trep jT'de 


the expres- 
Sm{nSuJnTrep) 


/"T/ 

V ^J-,multi) 


< 


_poo / / Tft_^Ezix,z)dz 


(A3) 


where M 


F°° = 


sillies UJnErep') 

2 COSh(fErep/Tdec^ COs(Jc^TiT^-gp)] 


(A4) 


We note that the expression (jA3l) can be used to find 
the deflection and time delay using the approach devel¬ 
oped for a single bunch. However, it requires detailed 
knowledge of E^ which can be obtained only with very 
accurate numerical analysis. To make reasonable estima¬ 
tions the following steps can be taken. One notes that 
if the driving bunches with transverse offset r from the 
mode electrical centre is acting upon trailing (by dis¬ 
tance z or time interval t = zfc (not r)) charges, it 
will provide a transverse kick via the excitation of the 
transverse dipole mode. Assuming a multipole represen¬ 
tation as discussed above (USD one can define dz =cdt 
and \7±[Ez{r,cj),z)] in (IA3I) as: 


M 

V±{r,t) = 2qr ^ A:„ 

m—1 


sin(a;mt)e =5" 


(A5) 
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where Wm, Qm are the frequency and the Q factor of mul¬ 
tipoles with index m (monopole modes are not consid¬ 
ered as they will not contribute to the kick). The ampli¬ 
tudes Kra of multipoles’ potential is also known as “kick 
factors” and are given by Km = c/4[i?/Q]™. Applying 
([2|) and taking into account only the first dipole mode 
m=l (i.e. ignoring high-order multipoles) Km can be 
presented as 


To study the damping requirements, which ensure the 
beam does not exceed the specified deflection, we under¬ 
take a single-mode analysis and hence we will drop the 
subscript m or, in case of (IA6I) . subscript 1. Following 
from ( 151171 ) . and applying expressions (IA21IA6I) observed 
to (fTTO) . we estimate radial offset between the beam 
trajectory and designed pass 


Ki 


V({x) 


dwir^C/i 


(A6) 


where r indicates the transverse deviation from the de¬ 
signed trajectories in accelerating (superscript j=l) and 
decelerating (superscript j=2) cavities, Um is the total 
energy of the mode, and oJm is the mode eigen-frequency. 


V^,max = 2IrKmF^luj (A7) 

Taking into account © with Ar = D 2 one finds that the 
expression for the maximum current: 


Ima. = {D2Wuj)/{2xKmF°^) (A8) 








